Journal of Computational Neuroscience

The phenomenon of splitting was originally observed in hamsters which, after prolonged exposure to constant light, exhibit two rest/wake cycles within a subjective day. Splitting is a consequence of the left and right suprachiasmatic nuclei (SCN) falling out of synchrony. While it is known that split activity is characterized by an antiphase relationship between the left and right SCN and between the core and shell within each hemisphere, the role of the commissural projections that connect the right and left SCN is not known. In the present study, we investigate the impact of the inter-hemispheric connections on the split and unsplit dynamics of a computational model of the bilateral SCN. Our model has 4 nodes corresponding to each right and left core and shell. We simulated our bilateral model under different lighting conditions and measured its period and the phase relationships among the 4 nodes. To further characterize the dynamics of the system, we performed a bifurcation analysis. We found that the bilateral model automatically splits unless entrained by bright light/dark cycles, or unless it has excitatory inter-hemispheric connections. This suggests that excitatory cross-connections may be important for freerunning behavior. We found that constant light of varying intensities transitions the model between split and unsplit activity only in very limited conditions, but the strength and polarity of the contralateral connections play a much greater role in this dynamical transition. These findings suggest that splitting may involve plasticity of the inter-hemispheric connections of the SCN.

Growing interest to describe the electrical behavior of glial cells, mainly astrocytes, in intact brain tissue poses more and more challenges to commonly accepted belief they only respond in a linear manner in uptake of the excess of extracellular potassium and maintenance of their network equipotentiality. Their highly conductive mutual interconnections via gap junction (GJ) connections introduce yet another class of nonlinear elements. As more studies report nonlinearities in membrane voltage Vm dependence of both, the membrane and junctional conductances, the need to formulate minimal dynamical models of their transient behavior is getting more acute. Since ODE models of coupled cells, even in simplest 1-d arrays, require simplified descriptions and small set of parameters, rare quantitative studies on glia makes the task even more difficult. This study attempts to qualify a self-coupled cell, or a glial cell coupled to fixed voltage as useful system for detecting the nature of instabilities and transitions coming from coupling. In a novel biophysical model of coupled astrocyte, we introduce nonlinear kinetics of deactivation for large junctional voltages for the first time. We found that N-shaped nonlinearities and corresponding fold structure in the vector field of isolated cell serves as a baseline on top of which coupling nonlinearities enrich the bifurcation picture. Numerical simulations of 1-d array of coupled astrocytes show that coupling increases the propensity of astrocytic Vm to bistability and front propagation. We believe that presented illustrations of possible effects of coupling nonlinearities will motivate neurobiologists to further explore their impact in disease. Significance statementTransient changes in membrane voltage of glial cells may produce significant transient voltage difference between directly coupled cells. Nonlinear steady-state conductance of their interconnection elements, the gap junctions, introduce nonlinear current profiles which are very difficult to measure and quantitate using the available methods due to marked permeability of the junctions and leakiness of glial membrane in general. We propose a minimal model of glial membrane extended with a self-coupled feedback loop, which under realistic simplifying assumptions could serve for qualitative analysis of the impact of coupling, on the stability of resting membrane voltage. Neuronal cells of the brain and spinal cord cannot exist and function without supportive and neuromodulatory functions of the diverse population of glial cells. This applies to virtually all physiological processes on cell level - from cell development, metabolic support, membrane signaling, slow molecular signal transduction, ion homeostasis, neurovascular coupling, myelination, to mention only a few, manifest neuro-glial interaction. Even though all glial cell types are interconnected, the most abundant ones, the astrocytes are massively interconnected by gap junctions to form ordered networks. Electrically, astrocytic networks display membrane voltage equipotentiality, which is considered system-wide resting state for given neuro-glial circuit or unit. With molecular and cellular substrates of glial connectivity being slowly elucidated, network science and dynamical modeling are slowly "invading" that area with many important issues left open. In this study using classical dynamical systems approaches we give indications how nonlinear intercellular coupling between astrocytes affects physiological resting state and its instabilities compared to isolated, uncoupled cell. We strongly believe the suggested minimal model could fill the gap in ODE modeling of neuro-glial circuits, within broadest scope of hypothesis-driven research in cell-level neuroscience.

The routing of information across the brains structural network is central to its wide range of functional capabilities. However, the mechanisms underlying information routing in complex brain networks, particularly between regions that do not share a direct anatomical connection, remain poorly understood. Neural mass models (NMMs), a computational modelling framework capable of capturing complex neural dynamics across scales, can potentially be used to study the dynamical and network bases of these vital polysynaptic routing processes. In this study, we investigate polysynaptic signalling in three widely used NMMs, obeying Ornstein-Uhlenbeck, Stuart-Landau, and Jansen-Rit dynamics, by tracking the propagation of a discrete, focal, high-amplitude perturbation across the underlying network. We find that polysynaptic propagation emerges in all tested NMMs when configured within dynamical regimes that effectively enhance the persistence of perturbations. We also find distinct parameter domains that maximise signal propagation to directly connected regions and to those separated from the source by at least two hops. Finally, we benchmark in silico stimulus propagation in the brain network against an empirical dataset of direct electrical stimulation trials, to explore the relative capabilities of the NMMs in capturing signal propagation to connected versus unconnected regions. This analysis highlights the significance of dynamical repertoire in capturing stimulus propagation outcomes. Overall, this study provides insights into how dynamical and network features shape signal propagation over complex brain networks.

The awake brain is known to display asynchronous (AS) states during periods of attention and arousal, but the responsiveness properties of such states remain unclear. Here, we investigate this question using computational models of spiking networks of excitatory and inhibitory neurons, mimicking recurrently-connected networks in layer 2/3 of the cerebral cortex. The networks can generate a continuum of AS states, but with different responsiveness characteristics. By using a mean-field model to infer the dynamic properties of the system, we find that there are two families of AS states, which we call "underdamped" (UD) and "overdamped" (OD). Responsiveness is maximised at the transition between OD and UD states, which identifies a "working point" that may present advantageous computational properties.

A growing number of studies indicate the possible involvement of astrocytes in triggering or modulating neurovascular coupling (NVC), i.e. the local dilation of blood vessels in the brain in response to neuronal activity. Astrocytes possess specialized subcellular compartments, named endfeet, that surround arterioles and capillaries, ideally positioned to mediate NVC. Various vasodilators have been shown to contribute to NVC, such as epoxyeicosatrienoic acid (EET), nitric oxide (NO), or prostaglandin E2 (PGE2), but the precise mechanisms underlying NVC and their variability remain to be fully elucidated. In particular, the involvement of astrocytes in this process is controversial. Recent translatome and proteomics data reveal that astrocytes and in particular endfeet are enriched in the proteins of the PGE2 pathway. However, how the latter could contribute to NVC remains to be characterized. Here, we develop a computational model of astrocyte-mediated NVC that recapitulates these findings and describes Ca2+ and PGE2 signaling in astrocytes, NO release by neurons, and arteriole diameter dynamics using ordinary differential equations. The model successfully reproduces the dynamics of arteriole diameter change during hyperemia from in vivo neocortical recordings in awake mice. Our simulations suggest that the astrocyte PGE2 pathway could be responsible for the late response of NVC at the arteriolar level. We further observe that PIP2-derived diacylglycerol plays a major role in driving arteriole diameter dynamics in our model, while phosphatidic acid-derived diacylglycerol, which is calcium-dependent, mainly acts as an amplifier of this response. Finally, a spatial implementation of the model using a simplified astrocyte geometry suggests that NVC is more efficient when synaptic stimulation occurs at the endfoot level rather than at other astrocytic compartments. Overall, this computational study suggests a partial role for astrocyte-mediated PGE2 release in NVC and points to astrocyte perivascular processes as sub-compartments that are ideally positioned and equipped to mediate NVC. Author summaryIn the brain, the local blood flow is regulated to meet neuronal energy demand by modulating the dilation of neighboring blood vessels. The mechanisms driving this process, known as neurovascular coupling (NVC), remain debated and are likely to differ depending on the physiological context. Recent evidence points to astrocytes, a cell type possessing specialized protrusions called "endfeet", that envelop the entire brain vascular tree. Contacts between synapses and endfeet have recently been reported, positioning the latter as ideal mediators of NVC. Here, we developed a computational model that simulates the signaling between neurons, astrocytes, and blood vessels. Our model successfully reproduces experimental recordings of blood vessels dilation in the brains of awake mice. Our simulations suggest that a specific signaling pathway in astrocytes, involving a molecule called prostaglandin E2, is a key driver of the late phase of NVC, occurring a few seconds after neuronal activity. Furthermore, our model indicates that the location of the stimulated synapses matters: signals sent to the astrocyte endfeet are particularly effective at controlling blood flow. This work helps clarify the active role of astrocytes in brain blood flow regulation, a process critical for healthy brain function.

Understanding how neuronal populations interact to encode and transform sensory information is a fundamental challenge in computational neuroscience. Most existing studies, however, study neural encoding, behavioral readout, and functional connectivity as disjoint problems. Two-photon calcium imaging enables simultaneous recording of large neuronal ensembles in vivo, driven by diverse stimuli and eliciting distinct behaviors. However, extracting directional functional connectivity metrics as well as encoding and readout properties of neurons from such data remains difficult due to indirect and noisy observations of spiking activity, slow temporal dynamics, and the latent interplay between external stimuli and endogenous neural processes. Here, we introduce a unified conceptual and operational modeling and inference framework for directly extracting functional Granger causal (GC) effects between neurons, from external stimuli to neurons, and from neurons to behavior, from two-photon imaging data, in the sense of Granger. Inspired by the intersection information framework, we also identify neurons that encode features of sensory stimuli that inform behavioral readout. The resulting GC networks together with the taxonomy of functional sensori-behavioral relevance, which we call G-taxonomy, provides a powerful statistical analysis framework, enabled by the integration of several techniques including state-space modeling and inference, variational inference, and point processes. We applied the proposed framework to simulated and experimentally-recorded two-photon imaging from the mouse auditory cortex (A1) during both passive listening and active tone discrimination. Our simulation studies reveal significant improvement of our proposed methodology over existing techniques. Analysis of experimental data from the mouse A1 identifies distinct groups of cells with diverse sensori-behavioral relevance, as well as changes in functional connectivity associated with correct vs. incorrect behavior. In summary, this work provides a principled and data-driven methodology for uncovering directional interactions among the neurons, sensory stimuli, and behavior, all within the same statistical framework, offering new insights into how distributed cortical populations transform sensory inputs into behaviorally relevant representations. Author SummaryThe brain processes sensory inputs through the coordinated activity of large networks of neurons and produces readouts that elicit behavior. Understanding how information flows and is processed through these networks is a central goal of neuroscience. In this study, we present a new computational framework that identifies directional interactions among neurons in an ensemble as well as from sensory stimuli to neurons and from neurons to behavior. Utilizing the Granger formalism to identify directional effects, as opposed to common correlational measures, our framework extracts said effects directly from two-photon calcium imaging data. We tested our proposed method on both simulated data and recordings from the auditory cortex of mice during passive listening and active tone discrimination tasks. Our method revealed diverse groups of neurons in the auditory cortex with distinct functional roles and relevance to sensori-behavioral integration. Our framework provides a new way to study the flow of information in the brain and can be broadly applied to uncover neural computations across sensory and cognitive systems.

Adolescence is an important developmental period during which there are diverse changes in the brain and behavior. Goal-directed behaviors and the component processes underlying those behaviors improve during adolescence, including working memory, response inhibition, and reinforcement learning. At the same time there is substantial pruning of excitatory connections in prefrontal cortex and ongoing myelination of axons. However, psychiatric disorders also become increasingly prevalent in late adolescence and early adulthood. In this study, we develop computational models that suggest a hypothesis for how the ongoing changes in the brain can give rise to the increased prevalence of psychiatric disorders. We show that both myelination and pruning during adolescence lead to attractor landscapes in which strongly encoded memories, driven by three-factor learning rules that modulate Hebbian plasticity, come to dominate the landscape of brain activity, at the expense of weakly encoded memories. Pruning and myelination lead to large, strong attractors which, if they are related to aversive emotions, can drive intrusive thoughts and compulsions in obsessive compulsive disorder, rumination in depression, and aversive memories in post-traumatic stress disorder. The link between pruning, myelination and the emergence of dominant attractors for emotionally salient memories is well supported by the models. The way these effects map onto forebrain circuits requires more work.

In previous studies by the author on binocular vision with the asymmetric eye (AE), which models a healthy human eye with misaligned optical components, the results were primarily presented in the Rodrigues vector (RV) framework and supported by simulations and 3D visualizations in GeoGebras dynamic geometry environment. In this paper, the novel geometric kinematics of the human eye, that is, the eye with misaligned optics, and simplified assumptions about the eye rotations (the eyes translational movements are disregarded), are developed within the framework of rigid-body rotations. The originality of the analysis lies in a precise geometric decomposition of a full rotation of the eyes posture into a torsion-free rotation (the geodesic part) and a torsional rotation (the non-geodesic extension of the geodesic part). This decomposition is extended to the corresponding decomposition of the angular velocity. A novel derivation of the eyes angular velocity from the RV formulation of the eye kinematics is proposed.

Electrical transmission is mediated by intercellular channels that cluster into structures known as gap junctions (GJ). In vertebrates, GJ channels are encoded by the gene family of connexin (Cx) proteins that assemble as hexamers, termed hemichannels, in the pre- and postsynaptic membranes, and that subsequently dock to form GJ channels. Auditory contacts on the fish Mauthner cells serve as model to study the properties and organization of vertebrate electrical synapses. Electrical transmission at these synapses is mediated by multiple co-existing GJs at which the presence of intercellular channels is regulated by a molecular scaffold. Zebrafish contain four homologs of the neuronal Cx36: Cx35.5 and Cx35.1 (gjd2a and b, respectively), and Cx34.1 and Cx34.7 (gjd1a and b). Cx mutations suggested that GJs are formed by heterotypic channels made of presynaptic Cx35.5 and postsynaptic Cx34.1. Using transgenic fish in which Cxs were tagged, we found that a second Cx, Cx34.7, is present together with Cx34.1 on the postsynaptic side at some but not all GJs at these terminals. When exogenously expressed, both Cx34.1 and Cx34.7 formed heterotypic functional channels with Cx35.5, each with substantially different voltage-dependent properties, indicating they can serve differential functions. However, we previously demonstrated that electrical transmission is lost in Cx34.1 but not Cx34.7 null mutants, suggesting that Cx34.7 cannot compensate for the loss of Cx34, despite the intrinsic ability of Cx34.1 and Cx34.7 to create functional channels. The findings reveal an unanticipated functional organization in the electrical synapse, where Cx34.1 is obligatory and Cx34.7 accessory, roles that appear to be defined by the postsynaptic molecular scaffold, with two postsynaptic Cxs possibly assembling under specific functional contexts. Thus, our results indicate that electrical synapses share an organizational motif with chemical synapses, akin to how they combine postsynaptic receptor types to modify synaptic function.

Understanding intracranial pressure (ICP) dynamics is essential for interpreting clinical infusion tests used in the diagnosis of cerebrospinal fluid circulation disorders. However, the complex coupling between vascular pulsations, cerebrospinal fluid flow, and intracranial compliance makes quantitative interpretation of these tests challenging. Here, I present a patient specific simulation framework based on an extended electrical analog model that reproduces intracranial pressure dynamics observed during clinical infusion tests. The model integrates physiological inputs including arterial blood pressure, heart rate, respiratory rhythm, and resistance to cerebrospinal fluid outflow derived from clinical data. Built upon the classical Ursino framework, the model incorporates several modifications enabling realistic representation of physiological pulsations and infusion test conditions. The resulting system functions as a hybrid electrical-numerical simulation model representing a simplified digital electrical twin of intracranial hydrodynamics. The model was validated using data from 21 clinical infusion tests performed in patients with suspected normal pressure hydrocephalus. Simulated intracranial pressure recordings were compared with clinical measurements using regression and residual analysis. The simulations demonstrated strong agreement with measured data, with a mean correlation coefficient of r = 0.95 (95% CI 0.94 - 0.96), mean residual values within -1.71 to +1.68 mmHg, and a mean root mean square error (RMSE) of 2.07 mmHg. These results demonstrate that the proposed model accurately reproduces the dynamic behavior of intracranial pressure observed during clinical infusion tests. The framework provides a physiologically grounded computational tool for studying patient specific intracranial dynamics and may support improved interpretation of infusion test results in clinical practice.

The learning dynamics of biological brains and artificial neural networks are of interest to both neuroscience and machine learning. A key difference between them is that neural networks are often trained from a randomly initialized state whereas each brain is the product of generations of evolutionary optimization, yielding innate structures that enable few-shot learning and inbuilt reflexes. Artificial neural networks, by contrast, require non-ethological quantities of training data to attain comparable performance. To investigate the effect of evolutionary optimization on the learning dynamics of neural networks, we combined algorithms simulating natural selection and online learning to produce a method for evolutionarily conditioning artificial neural networks, and applied it to both reinforcement and supervised learning contexts. We found the evolutionary conditioning algorithm, by itself, performs comparably to an unoptimized baseline. However, evolutionarily conditioned networks show signs of unique and latent learning dynamics, and can be rapidly fine-tuned to optimal performance. These results suggest evolution constitutes an inductive bias that tunes neural systems to enable rapid learning.

IntroductionAlthough neural activity is organized across temporal and spatial scales, the principles that determine the accuracy and fidelity of neural information representation across scales remain unclear. In particular, while recent empirical results have reported mesoscopic optimality in neural decoding, no theoretical accounts exist that explain when and why such intermediate scales emerge as optimal. Here, we develop an analytical framework to study the optimal temporal scale of information representation and its dependence on the dynamic structure of signal and noise in neural data. Materials and MethodsWe formulate a multiscale theoretical model in which neural population activity is represented by temporally encoded trial vectors at microscopic, mesoscopic, and macroscopic resolutions. Neural responses are modeled as class-dependent mean activations (signal) corrupted by temporally correlated noise, and decay rates of correlations in both signal and noise are varied parametrically. Representational quality at each scale is quantified using the sensitivity index (d-prime) for decoding condition from neural activity. ResultsWe derive closed-form expressions for the sensitivity index at each temporal scale. These expressions reveal the key roles of signal and noise correlations as the main determinants of condition decodability at all scales. Comparing expressions under various combinations of signal and noise correlations reveals two regimes. First, when signal and noise correlations are absent or persistent over time, the optimal resolution falls at one of two extremes: macroscale (resp. microscale) if signal correlations are stronger (resp. weaker) than noise correlations. In contrast, when both signal and noise correlations decay with temporal separation, temporal integration produces a nontrivial trade-off: moderate integration improves decodability by suppressing noise while preserving coherent signal, whereas excessive integration degrades signal and decodability. Therefore, only in the latter regime, mesoscopic representations emerge as optimal across a broad range of biologically plausible parameters. DiscussionThis work provides a theoretical explanation for how the optimal temporal scale of neural information representation depends on the interplay between signal and noise correlations and their temporal decay. Broadly, the framework establishes temporal integration as a principled mechanism linking multiscale neural dynamics to information representation and offers testable predictions across recording modalities and neural systems.

Rhythm is a key building block of human music, speech and numerous other human activities. Understanding the computational substrates of rhythm perception requires models that bridge algorithmic function with biological implementation. We propose a physiologically grounded spiking neural network (SNN) framework to investigate the emergent representation and interpretation of auditory rhythms. Utilizing a recurrent SNN architecture trained on an auditory entrainment task, we characterize the networks latent dynamics through the analysis of firing rates and membrane potential fluctuations. Our results demonstrate that simulated neural populations exhibit phase-locking to the stimulus beat, with endogenous oscillations driven by rhythmic input. We further show that anticipatory dynamics--characterized by pre-stimulus depolarization--emerge naturally from the networks synaptic plasticity and temporal integration properties, rather than from explicitly defined oscillators. By treating network layers as functional analogs of cortical populations, this framework allows for the application of spectral and information-theoretic analyses typical of empirical electrophysiology. More in general, this approach establishes SNNs as robust exploratory tools for uncovering how predictive coding and rhythmic entrainment arise from the inherent constraints of biological neural computation.

The Poisson-Nernst-Planck (PNP) system is an accurate model of electrodiffusion of ionic species. It is commonly used in situations where nanoscale resolution is required, for instance close to ion channels in the membranes of biological cells. The inherent stiffness of the equations has made them challenging to solve and has limited the applicability of the system. In particular, the time step required for stable solutions has typically needed to be very short (nanoseconds), which makes simulations on the time scale of an action potential (milliseconds) difficult. Recently, it has been observed that avoiding operator splitting and instead solving the concentration equations and the electrostatic equation in a coupled manner relaxes the time-step limitation considerably. However, no theoretical explanation of this observation has been provided. Here, we aim to explain why the coupled scheme allows much larger time steps. We illustrate the mechanism by considering special cases that define necessary, but not sufficient, conditions for stability. We also show that these conditions remain relevant for the fully coupled PNP model in 3D.

We present a free and open-source, semi-automated, topologically robust pipeline for fitting cable models to 3D surface mesh morphology data of neuronal membranes, particularly suited to structures with complex shapes and topological holes. The motivation for this work is the discovery of morphologically complex neural spines on the auditory space-specific neurons of the barn owl (Tyto alba, Tyto furcata), dubbed "toric spines", notable for their high curvature, branching density, and holes/loops. Multicompartmental simulation software requires morphology to be represented as cable models (e.g., SWC format), yet existing software tools for fitting cable models to complex 3D surface meshes have not produced satisfactory results for toric spines, and loops are generally unsupported. We present the Mesh and Skeleton Cable Fitting (MASCAF) pipeline and software, which fits a cable model (e.g., SWC format) to a surface mesh using mean-curvature flow skeletonization. In this paper, we demonstrate how MASCAF is applied to fit cable models, how loops can be reconstructed in simulations with the Arbor and NEURON simulation software, and how the results can be validated using geometry and simulator-based methods. While non-tree morphologies such as toric spines are neuroanatomically special, our software pipeline provides a cable-model fitting approach for surface mesh data that is topologically robust, deterministic, open-source, and applicable to general morphologies, thereby closing a crucial gap between neuronal imaging and high-resolution simulation.

A motoneuron pool is often regarded as a rigid controller because the largely shared synaptic input across motoneurons leads to strongly correlated activity. However, brief deviations from this correlated behavior have been observed even in some constrained tasks, raising the question of whether these results reflect limitations of the rigid view of motoneuron pool control. Here we show that they do not. We developed a biophysical model of a motoneuron pool receiving shared excitatory and inhibitory synaptic inputs that also included the motoneuron-specific effects of neuromodulation; model parameters were tuned based on large-scale motoneuron recordings in humans. Simulations showed that the intrinsic differences in how motoneurons respond to neuromodulation are both necessary and sufficient to transiently decorrelate pairs of motoneurons receiving a shared synaptic input. Crucially, such transient decorrelation is only observed when motoneurons have different sensitivity to neuromodulation, consistent with experimental observations during volitional control in humans. Our model also explains how participants can improve their ability to transiently decorrelate the activity of motoneurons innervating the same muscle by leveraging refined behavioral strategies that exploit the differential response of motoneurons to neuromodulation, rather than through physiological changes. These results identify that heterogeneous sensitivity to neuromodulation enables brief flexibility in the otherwise rigid control of motoneurons enforced by a shared synaptic input, and show how practice allows participants to exploit latent flexibility within otherwise rigid constraints.

Vestibular neurons are core elements of the pathways involved in vestibulo-motor functions, such as vestibulo-spinal and vestibulo-ocular reflexes. To meet behavioral needs, electrophysiological neuronal properties are adequately adapted to the sensory-motor computation sustaining these distinct vestibular reflexes. During frog metamorphosis, there is a complete reorganization of the posturo-locomotor system while the oculomotor system remains minimally changed, probably associated to so far unknown changes in vestibular neuronal properties. We used this unique model to investigate the central developmental mechanisms underlying such a reconfiguration of vestibular-associated behaviors. Central vestibular neurons exhibit two types of electrophysiological phenotypes: tonic neurons with a continuous discharge and phasic neurons with a transitory discharge mainly due to the activation of Kv1.1 channel. Electrophysiological recordings and Kv1.1 immunolabeling of vestibulospinal (VS) and vestibulo-ocular (VO) neurons at both larval and juvenile stages revealed that the majority of VS neurons exhibited a tonic discharge in larvae but a phasic discharge in juvenile, while VO neurons remained mainly tonic throughout development. Changes in phasic and tonic neurons proportions in VS population are partly explained by neurogenesis. But we provide evidences that an electrophysiological phenotype switch is a concomitant developmental mechanism participating in the maturation of these central vestibular neurons. All together our results showed that the maturation process in central vestibular neuronal groups is highly related to the metamorphosis-induced remodeling of vestibulo-motor functions they are involved in, with the ultimate purpose of ensuring an adequate adaptation of neuronal elements properties to the developmental changes of behavioral constrains.

BackgroundThe relevance of spontaneous activity has been unlocked thanks to recent large scale recordings that revealed, via Shared Variance Component Analysis (SVCA), the high-dimensional nature of the ongoing activity. A fundamental problem is how the dimension modifies when more neurons are included in the analysis. Contradictory results have been reported on this subject based on SVCA and Principal Component Analysis (PCA). New MethodWe investigate pro et contra of SVCA and PCA for the identification of reliable responses encoding underlying state variables. We focus on common features of the spectra of the reliable variances (RVs) and on their dimensionality. The analysis is demonstrated on previously published Ca2+ data from the visual and the dorsal cortex in head fixed mice during spontaneous behavior. ResultsRVs grow proportionally to the number N of neurons and show a power-law decay k- with the k-th SVC dimension over a range bounded by a maximal dimension kc, initially diverging as N 1/ and then saturating at sufficiently large N. The reliable dimensionality, estimated with different methodologies, also shows a clear saturation to an asymptotic value for large N. Furthermore, its value decreases when becomes larger, as demonstrated by employing experimental data as well as theoretical predictions. ConclusionWe have shown that SVCA is an extremely effective tool to extract reliable features from the neural signals, and that the exponent represents a biomarker able to reveal the level of correlation of the neurons as well as the dimensionality of the reliable space. HighlightsO_LIAdvantages and drawbacks of Shared Variance Component Analysis to extract reliable signals from neural data C_LIO_LIComparison of different methods to estimate reliable neural dimensionality associated to spontaneous activity C_LIO_LIAnalytical expressions of embedding dimensionality for power-law decaying reliable variances C_LIO_LIBounded growth of the dimensionality with the number of neurons C_LI

BackgroundThis study examines a competition-based model (C-model) designed to capture the temporal dynamics of successive brain microstates derived from electroencephalography (EEG) recordings during eyes-open conditions. The analyzed data were obtained from a public repository comprising microstate sequences from 60 sessions of a single subject [1]. When applied to microstate dynamics, the C-model posits a stochastic competition among neural circuits underlying the expression of individual microstates. MethodsThe model is formulated at a conceptual level (computational level in Marrs framework) and employs a geometric distribution to account for the long right tail of microstate duration distributions, interpreted as the probability of "failure" of the currently active microstate to persist. To account for the short-lived left tail, the model incorporates a transient increase in the stability of the currently active network, or equivalently, a temporary decrease in the activation probability of competing microstates (refractory period). ResultsThe model provides a good fit to the microstate duration distributions across all 60 sessions. One third of sessions showed microstate identity sequential dependency with respect to the previous microstates. DiscussionThese results suggest that the C-model captures key aspects of microstate temporal structure. Moreover, because microstate probabilities can be modulated by psychophysiological conditions--including the influence of previously active networks--the model may serve as a building block for more comprehensive neurobiological frameworks of neural and behavioral dynamics. In such frameworks, microstate sequences could emerge from structured competition and flow among neural networks supporting microstate expression.

Neuronal computation depends on the balance between excitation and inhibition, yet how this balance is implemented across the dendritic tree remains unclear. Classical views predict that inhibition should be most effective near the soma or along the path from excitation to output, but many interneuron subtypes preferentially target remote dendritic compartments. This apparent paradox is sharpened by active dendrites, where local NMDA spikes, calcium plateaus and backpropagating action potentials can make distal branches powerful contributors to somatic firing. Here we develop an analytical framework that extracts general principles of inhibition from biophysically detailed multi-compartment simulations. By reformulating the implicit voltage update of detailed neuron models as a matrix recursion, we derive exact voltage sensitivities to inhibitory synaptic perturbations. This leads to a unified {Phi}-a law: the somatic impact of inhibition factorizes into a global dendritic susceptibility term and a local synaptic perturbation term. Using this law to map inhibitory leverage and identify optimal inhibitory interventions, we show that active dendritic excitation can shift inhibitory hot zones from perisomatic regions toward distal or intermediate compartments. Across neocortical, hippocampal and striatal neuron models, the same response law explains convergent inhibitory strategies despite distinct cellular mechanisms. Our framework turns detailed numerical simulation into analytical theory, providing a general principle for how diverse dendritic inhibition controls active neurons.